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Abstract 



The LHC program will include the identification of events with single prompt high-/cj_ photons as 
probes of new physics. We show that this channel is uniquely suited to search for experimental 

J*. evidence of TeV-scale open string theory. At the parton level, we analyze single photon production 

in gluon fusion, gg — > jg, with open string states propagating in intermediate channels. If the 
photon mixes with the gauge boson of the baryon number, which is a common feature of D-brane 
quivers, the amplitude appears already at the string disk level. It is completely determined by 

j> . the mixing parameter (which is actually determined in the minimal theory) - and it is otherwise 

model- (compactification-) independent. We discuss the string signal cross sections as well as the 



QCD background. The present analysis takes into account the recently obtained decay widths of 
first Regge recurrences, which are necessary for more precise determination of these cross sections 
in the resonant region. A vital part of the background discussion concerns the minimization of 



misidentified n s emerging from high-pj_ jets. We show that even for relatively small mixing, 



100 fb _1 of LHC data could probe deviations from standard model physics associated with TeV- 
scale strings at a 5<r significance, for M s t r i ng as large as 2.3 TeV. It is also likely that resonant 
k> i bumps could be observed with approximately the same signal-to-noise ratio. 



I. GENERAL IDEA 

The CERN's Large Hadron Collider (LHC) is the greatest basic science endeavor in 
history. Spectacular physics results are expected to follow in short order once it turns on 
this year. LHC will push nucleon-nucleon center-of-mass energies up to \/s = 14 TeV for 
pp collisions and yfs = 5.5 TeV for Pb-Pb collisions. The ATLAS and CMS detectors 
will observe the highest-energy particle collisions produced by the accelerator, whereas the 
ALICE detector will observe the very messy debris of heavy ion collisions. The LHC will 
probe deeply into the sub-fermi distances, committing to careful searches for new particles 
and interactions at the TeV scale. 

At the time of its formulation and for years thereafter, Superstring Theory was regarded 
as a unifying framework for Planck-scale quantum gravity and TeV-scale Standard Model 
(SM) physics. Important advances were fueled by the realization of the vital role played 
by D-branes |l| in connecting string theory to phenomenology (2j. This has permitted the 
formulation of string theories with compositeness setting in at TeV scales |3j and large 
extra dimensions. There are two paramount phenomenological consequences for TeV scale 
D-brane string physics: the emergence of Regge recurrences at parton collision energies 
Vs~ string scale = M s ; and the presence of one or more additional U(l) gauge symmetries, 
beyond the U(l)y of the SM. The latter follows from the property that the gauge group for 
open strings terminating on a stack of N identical D-branes is U(N) rather than SU(N) for 
N > 2. (For N = 2 the gauge group can be Sp(l) rather than U(2).) In this paper we exploit 
both these properties in order to obtain a "new physics" signal at LHC which, if traced to 
low scale string theory, could with 100 fb _1 of data probe deviations from SM physics at a 
5cr significance for M s as large as 2.3 TeV. A short version highlighting the salient results 
of our analysis has been issued as a companion Letter [4|. The present analysis, however, 
takes into account the recently obtained decay widths of first Regge recurrences [5j, which 
are necessary for more precise determination of cross sections in the resonant region. 

To develop our program in the simplest way, we will work within the construct of a 
minimal model in which we consider scattering processes which take place on the (color) U(3) 
stack of D-branes. In the bosonic sector, the open strings terminating on this stack contain, 
in addition to the SU(3) octet of gluons, an extra U(l) boson (C^, in the notation of J6|]), 
most simply the manifestation of a gauged baryon number symmetry. The U(l)y boson Y^, 
which gauges the usual electroweak hypercharge symmetry, is a linear combination of C^, the 
£7(1) boson B^ terminating on a separate U(l) brane, and perhaps a third additional U(l) 
(say W/j) sharing a U(2) brane to which are also a terminus for the SU(2)l electroweak gauge 
bosons W°. Thus, critically for our purposes, the photon A^, which is a linear combination 
of Y^ and W^ will participate with the gluon octet in (string) tree level scattering processes 
on the color brane, processes which in the SM occur only at one-loop level. Such a mixing 
between hypercharge and baryon number is a generic property of D-brane quivers, see e.g. 
Refs.pl, LZj, [sj] . The vector boson Z' , orthogonal to the hypercharge, must grow a mass M z > 
in order to avoid long range forces between baryons other than gravity and Coulomb forces. 

The process we consider (at the parton level) is gg — > gj, where g is an SU(3) gluon 
and 7 is the photon. As explicitly calculated below, this will occur at string disk (tree) 
level, and will be manifest at LHC as a non-SM contribution to pp — > 7 + jet. A very im- 
portant property of string disk amplitudes is that they are completely model- independent; 
thus the results presented below are robust, because they hold for arbitrary compactifica- 
tions of superstring theory from ten to four dimensions, including those that break super- 




FIG. 1: Open string disk diagram for gg — > (77 scattering. The dots represent vertex insertions of 
the gauge bosons on the boundary of the world sheet. 

symmetry. The SM background for this signal originates in the parton tree level processes 
gq — > jq, gq — > jq, and qq —* jg. Of course, the SM processes will also receive stringy 
corrections which should be added to the pure bosonic contribution as part of the sig- 
nal [9|, |lOj, lll|, ll2j|- We leave this evaluation to a subsequent publication [13J; thus, the 



contribution from the bosonic process calculated here is to be regarded as a lower bound to 
the stringy signal. It should also be stated that, in what follows, we do not include effects 
of Kaluza-Klein recurrences due to compactification. We assume that all such effects are in 
the gravitational sector, and hence occur at higher order in string coupling [9]. 

The plan of the paper is as follows. In the next section we outline the calculation of 
the string amplitude for the process pictured in Fig. [T] and show that there is no amplitude 
containing the zero mass poles of the SM. In Sec. IIHI we first calculate cross sections for 
gluon fusion in the resonance region as well as QCD background, for a simple fc±, m in cut 
on the transverse momentum of the photon. A vital part of the background discussion 
concerns the minimization of misidentified 7r°'s from high-p^ jets. After that we calculate 
the signal-to-noise ratio and show that a significant deviation from SM can be obtained for 
fc±,min > 300 GeV for M s as large a 2.3 TeV. In Sec. [IV] we delineate the search for resonant 
structure in the data. Our conclusions are collected in Sec. [V] The Appendices contain 
some additional formulae referred to in the main text. 



II. THE STRING AMPLITUDE 

The most direct way to compute the amplitude for the scattering of four gauge bosons 
is to consider the case of polarized particles because all non- vanishing contributions can be 
then generated from a single, maximally helicity violating (MHV), amplitude - the so-called 
partial MHV amplitude [14|. Assume that two vector bosons, with the momenta k\ and &2, 
in the U(N) gauge group states corresponding to the generators T ai and T a2 (here in the 
fundamental representation), carry negative helicities while the other two, with the momenta 
&3 and &4 and gauge group states T° 3 and T a4 , respectively, carry positive helicities. (All 
momenta are incoming.) Then the partial amplitude for such an MHV configuration is given 



by [15|,llfl 



A(l-,2~,3 + ,4+) = 4c/ 2 Tr(T ai T a2 T a3 T a4 ) / w \ / - - V(h, k 2 ,k 3 , fe 4 ) , (1) 

v ' ' ' ; y v ; (12) (23) (34) (41) v ' ' 3 ' ; ' v ; 

where g is the U(N) coupling constant, (ij) are the standard spinor products written in the 
notation of Refs. [121, ll8(, and the Veneziano formfactor, 



V{k u k 2 , k 3 , k 4 ) = V(s, t, u) = ^ r fj^ ~ , ( 2 ) 

is the function of Mandelstam variables, here normalized in the string units: 

2k 1 k 2 2k\k?, 2k\k± . 

s = ^p t = ^p M = ^T : * + * + - = °- ( 3 ) 

(For simplicity we drop carets for the parton subprocess.) Its low-energy expansion reads 

7T 2 

V(s,t,u) « 1- — su-((3)stu+... (4) 

6 

We first consider the amplitude involving three SU(N) gluons gi, g 2 , g% and one U(l) 
gauge boson 74 associated to the same U(N) quiver: 

rpai = T a ^ T a 2 = rpb ^ rpas = rpc ^ rpa^ = QJ ^ (5) 

where I is the NxN identity matrix and Q is the U(l) charge of the fundamental represen- 
tation. The U(N) generators are normalized according to 

Tr(T a T 6 ) = -5 ab . (6) 



Then the color factor 



Tr(T ai T a2 T a3 T a4 ) = Q(d abc + -f abc ) , (7) 

where the totally symmetric symbol d abc is the symmetrized trace while f abc is the totally 
antisymmetric structure constant. 



The full MHV amplitude can be obtained [l_5|, [16] by summing the partial amplitudes 



with the indices permuted in the following way: 

am - - + +1 / 2„n>4 V- Tr ( 7--r"-7--r") V(k u ,k 2 „,k 3 „k t ) 

where the sum runs over all 6 permutations a of {1,2,3} and i a = oil). Note that in the 
effective field theory of gauge bosons there are no Yang-Mills interactions that could generate 
this scattering process at the tree level. Indeed, V = 1 at the leading order of Eq.fTjJ and 
the amplitude vanishes due to the following identity: 

1 1 1 

+ /oo\/oiWi„W,o\ + /oiwowo„w„o\ = " ( 9 ) 



(12) (23) (34) (41) (23) (31) (14) (42) (31) (12) (24) (43) 



Similarly, the antisymmetric part of the color factor ((Tj) cancels out in the full amplitude 
(JSJ) . As a result, one obtains: 

M(g7,92,9t,lt) = 8Q d abc g 2 (12) 4 ( -J&^L— + __^£iM^| , (10) 
wi,y2,yz,uj v y\ I ^(12) (23) (34) (41) (12) (24) (13) (34) ^ ' v ; 

where 

, (s , t ,„) = r(i- H )(£|^-Q^|). (ii) 

All non-vanishing amplitudes can be obtained in a similar way. In particular, 

M(g7,gt,g^,it) = 8Qd a V(13) 4 ( , ^fi^r, : + , w^w^w r l , (12) 
Wi,y2,y3,uj ^ y\ 1 ^(13) (24) (14) (23) (13) (24) (12) (34) ^ ' v ; 

and the remaining ones can be obtained either by appropriate permutations or by complex 
conjugation. 

In order to obtain the cross section for the (unpolarized) partonic subprocess gg — > #7, 
we take the squared moduli of individual amplitudes, sum over final polarizations and colors, 
and average over initial polarizations and colors. As an example, the modulus square of the 
amplitude ([8]) is: 



\M(g^g 2 ,gt,"ft)\ 2 = 64Q 2 d abc d abc g 



sfj,(s,t,u) Sfi(s,U,t) 



u 



(13) 



Taking into account all 4(iV 2 — l) 2 possible initial polarization/color configurations and the 
formula [l9( 

a,b,c 

we obtain the average squared amplitude 



\M{gg^ gi )\ 2 = g A Q 2 C{N) 
where 



sfi(s,t,u) sfi(s,u,t) 



u t 



+ {s <-► t) + (s <-»■ u) } , (15) 



C ( N) = 2 f~\ (16) 

V ; N(N 2 -l) V ; 

The two most interesting energy regimes of gg — > gj scattering are far below the string 
mass scale M s and near the threshold for the production of massive string excitations. At 
low energies, Eq. ffT5|) becomes 

\M(gg -> 91 )\ 2 « g 4 Q 2 C(N)^(s 4 + 1 4 + u 4 ) (s, t, u « 1) . (17) 

The absence of massless poles, at s = etc., translated into the terms of effective field theory, 
confirms that there are no exchanges of massless particles contributing to this process. On 
the other hand, near the string threshold s ~ M 2 (where we now restore the string scale) 

\M(gg - gi )\ 2 « ^Q 2 C(N)^±^±^ 2 (s « M 2 ). (18) 



The singularity at s = M 2 needs softening to a Breit-Wigner form, reflecting the finite 
decay widths of resonances propagating in the s channel. Due to averaging over initial 
polarizations, Eq. (fl8|) contains additively contributions from both spin J = and spin 
J = 2 gluonic Regge recurrences (G* in the notation of Ref.[5J]), created by the incident 
gluons in the helicity configurations (±±) and (±=f), respectively. The M| term in Eq. (TTSJ) 
originates from J = 0, and the t 4 + w 4 piece reflects J = 2 activity. Since the resonance 
widths are spin-dependent [5(: 



r J=0 = -a s M s « 75 (M s /TeV) GeV . 
T J=2 = ^a s M s « 45 (M s /TeV) GeV 



(19) 



the pole term ffT8|) should be smeared as 



|A*(00->07)I : 



4 9 4 g 2 c(#) 

M 4 



M„ 8 



t 4 + u 4 



(s - M 2 ) 2 + (r J=0 M,) 2 (s - M 2 ) 2 + (r J=2 M s ) 2 



• (20) 



a. 



In what follows we will take N = 3 and set g equal to the QCD coupling constant, 
= (g 2 /47r) ~ 0.1. Before proceeding with numerical calculation, we need to make precise 
the value of Q. If we were considering the process gg — > C°g, where C° is the U(l) gauge field 
tied to the U(3) brane, then Q = a/1/6 due to the normalization condition ([6]). However, for 
gg — > 'jg there are two additional projections: from C^ to the hypercharge boson Y^, giving 
a mixing factor k; and from Y^ onto a photon, providing an additional factor cos 9w {&w 
Weinberg angle). The C°—Y mixing coefficient is model dependent: in the minimal model 
it is quite small, around k ~ 0.12 for couplings evaluated at the Z mass, which is modestly 
enhanced to k ~ 0.14 as a result of RG running of the couplings up to 2.5 TeV. It should 
be noted that in models [7|, [8J possessing an additional U(l) which partners SU(2) L on a 
U(2) brane, the various assignment of the charges can result in values of k which can differ 
considerably from 0.12. In what follows, we take as a fiducial value k 2 = 0.02. Thus, if ( 120]) 
is to describe gg — > jg, we modify our definition of Q given in Eq. (jSJ) to accommodate the 
additional mixings, and obtain 



Q 2 = 1 K 2 cos 2 fl w ~ 2.55 x 10~ 3 (k 2 /0.02) 



(21) 



In the remainder of the paper, we explore potential searches for Regge excitations of funda- 
mental strings at LHC. 

III. ISOLATED HARD PHOTONS 



In order to assess the possibility of discovery of signal above QCD background, we adopt 
the kind of signal introduced in [20] to study detection of TeV-scale black holes at the 
LHC, namel y a high-fc^ isolated 7 or Z. Thus, armed with parton distribution functions 
(CTEQ6D) [21] we have calculated integrated cross sections a(pp — > 7 + jet)|fc ± ( 7 )>fc Xmin 
for both the background QCD processes (see Appendix I) and for gg — > jg, for an array 
of values for the string scale M s (see Appendix II). Our results are shown in Fig. [3 It is 
evident that the background is significantly reduced for large k ±rain . At very large values 
of k± min , however, event rates become problematic. In Fig. [3] we show the string cross 
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FIG. 2: Behavior of the QCD cross section for pp — > 7 + jet (dot-dashed line) as a function of 
fc_Lmin- The string cross section overlying the QCD background is also shown as a solid line, for 
M. = 1 TeV. 



section and number of events (before cuts) in a 100 fb _1 run at LHC, for k±^ min = 300 GeV, 
as a function of the string scale M s . Next, we explore the LHC discovery potential by 
computing the signal-to- noise ratio ( signal /y/SM background = S/N). For a 300 GeV cut 
in the transverse momentum, the QCD cross section (shown in Fig. [2]) is about 8 x 10 3 fb, 
yielding (for 100 fb -1 ) \/SM background ~ 895. A point worth noting at this juncture: to 
minimize misidentification with a high-fc^ vr°, isolation cuts must be imposed on the photon, 
and to trigger on the desired channel, the hadronic jet must be identified [22J. We will leave 
the exact nature of these cuts for the experimental groups, and present results for a generous 
range of direct photon reconstruction efficiency. To do so, we define the parameter 



background due to misidentified ir° after isolation cuts 

(3 = h 1 

QCD background from direct photon production 



(22) 



Therefore, the noise is increased by a factor of \f]3, over the direct photon QCD contribution. 
Our significant results are encapsuled in Fig. HJ where we show the discovery reaches of the 
LHC for different integrated luminosities and k 2 = 0.02. A detailed study of the CMS 
potential for isolation of prompt-7's has been recently carried out [23j, using GEANT4 
simulations of 7 + jet events generated with Pythia. This analysis (which also includes 7's 



X3 



C 

.2 

o 

cu 

CO 

CO 
CO 

o 
O 



ov 



D ^ 



0°E 



Z E 



E 



1 E 



10 



\ 



\ 



\ 



k x , min = 300 GeV 
_._._ PbPb 
PP 



\ 



1.5 



2.5 



3.5 

M s (Tc 



>V 



4— 

o 

O 



to 



> 



o 

C 
X, 



FIG. 3: Cross section for gluon fusion into 7 + J e t|fc ± ( 7 )>3oo GeV an d expected number of events, 
for 100 fb~ and varying string scale. 

produced in the decays of 77, K®, u>°, and bremsstrahlung photons emerging from high-p^ 
jets) suggests (3 ~ 2. Of course, considerations of detector efficiency further reduce the S/N 
ratio by an additional factor e, where 1 < e ^C \J~ft. We conclude that discovery at the 
LHC would be possible for M s as large as 2. 3 Te V. The dependence of the discovery reach 
with the C° — Y mixing coefficient k has been extensively discussed in the accompanying 
Letter [j. 

We now briefly explore the potential of ALICE to search for low mass string excita- 



tions [24j]. With this motivation, we extend our analysis to include heavy ions collisions. In 



the spirit of Ref. [25j we consider the unshadowed parton distribution functions, i.e. 



R, 



i/A 



. X 



f i/A (x,Q) 
Afi(x,Q) 



(23) 



where fi/A and fi are the parton distribution functions inside a free nucleus of mass A 
and free nucleon, respectively. For M s > 1 TeV, this approximation holds because LHC 
Pb-Pb collisions probe the minimum value of parton momentum at x m ; n ~ M^/s ~ 0.033, 
where there are no shadowing effects. A comparison of the string cross section for gluon 
fusion into 7 + jet|fc ± ( 7 )>3oo GeV for pp and Pb-Pb collisions is shown in Fig. [31 However, the 
larger aggregate of partons also increase the SM background; namely, for fcj_ )m i n > 300 GeV, 
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FIG. 4: Contours of 5a discovery in the (M s , detector efficiency) plane for different integrated 
luminosities and k? = 0.02. 



cpb-Pb^7X ~ 2.8 x 10 7 fb. This greatly decreases the sensitivity to D-brane models, which 
would require a Pb-Pb integrated luminosity of a few hundred pb -1 . This is substantially 
larger than the present day estimate f2( 



IV. BUMP-HUNTING 



The discovery trigger described in the previous section, the observation of isolated pho- 
tons at large transverse momentum, serves very well as a signature of new physics. Given the 
particular nature of the process we are considering, the production of a TeV-scale resonance 
and its subsequent 2-body decay, signatures in addition to large k± photons are available. 
Most apparently, one would hope that the resonance would be visible in data binned accord- 
ing to the invariant mass M of the photon + jet, setting cuts on photon and jet rapidities, 
Vi, 2/2 < 2/max, respectively. With the definitions Y = \{yi + 2/2) and y = \(y\ — 2/2), the 
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FIG. 5: da/dM (units of fb/GeV) vs. M (TeV) is plotted for the case of SM QCD background 
(dashed) and (first resonance) string signal + background (solid). 



cross section per interval of M for pp — > 7 + jet + X is given by [2 

l '°- = Mr Yl I d\ ./,-(./■„. XI) i ;(,,;,. M) ! 

-(Vmax+y) 



dM 



MtJ2 dY f&a, M) frfa, M) dy — 

ijk L^-Vmax ; - 



d< 



ij—fyk 



cosh y 



Jmax Py max * A/x 

dY fi(x a ,M) fj(x b ,M) / dy — 

J — (j/max— V) 



dt 



1 



ij—+jk 



cosh y 



(24) 



where i, j, fc are different partons, r = M 2 /s, x a = \fre Y , and x^ = yfre Y . The kinematics 
of the scattering provides the relation 

fc = ;Ar (25) 

2 coshy 

which, when combined with the standard cut k± > &;_i_ jm i n , imposes a lower bound on y to be 
implemented in the limits of integration. (For details see Appendix III.) The Y integration 
range in Eq. (1241) . Y max = min{ln(l/ v /r), y max }, comes from requiring x a , x b < 1 together 
with the rapidity cuts \y\\, l^l < 2.4. Finally, the Mandelstam invariants occurring in the 
cross section are given by 's = M 2 , t = —\M 2 e~ y / coshy, and u = —\M 2 e +y / coshy. 
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In Fig. [5] we show several representative plots of this cross section for different values of 
M s . Standard bump- hunting methods, such as calculating cumulative cross sections 

r°° rlrr 

a(M ^L,iM dM w 

and searching for regions with significant deviations from the QCD background, may allow 
to find an interval of M suspected of containing a bump. With the establishment of such 
a region, one may calculate a signal-to-noise ratio, with the signal rate estimated in the 
invariant mass window [M S — 2T, M S +2T]. This estimate of signal-to-noise would be roughly 
the same as that obtained through the inclusive cut k± > 300 GeV. This follows from the 
relation ( |25|) : for M in the range of M s > 2 and for the significant contributing regions of y, 



the resulting k± cut in Eq. (|25|) does not differ significantly from the 300 GeV 28] . Should 
bumps be found, the D-brane model can be further differentiated from other TeV-scale 
resonant processes by the details of the angular distributions inherent in Eq. (I201) . 

V. CONCLUDING DISCUSSION 

In this work we have described how to search for the effects of Regge excitations of 
fundamental strings at LHC collisions. The underlying parton process for the excitation of 
the string resonance is dominantly the single photon production in gluon fusion, gg — > jg, 
with open string states propagating in intermediate channels. If the photon mixes with 
the gauge boson of the baryon number, which is a common feature of D-brane quivers, the 
amplitude appears already at the string disk level. It is completely determined by the mixing 
parameter - and it is otherwise model- (compactification-) independent. We have shown that 
even for relatively small mixing, 100 fb _1 of LHC data (in the pp — > 7+ jet channel) could 
probe deviations from SM physics at a 5<r significance, for M s as large as 2.3 TeV. We note 
that such a numerical value for the discovery reach is lower than the estimate presented 
in the accompanying Letter, M s ~ 3.3 TeV [4j. The present analysis contains a refined 
treatment of the resonance region, including the recently computed decay widths of both 
spin J = and spin J = 2 Regge recurrence of the gluon octet [5(. The discovery reach 
is lower because these resonances are slightly wider than naively expected while the signal 
cross sections are very sensitive to the width values. 

In closing we discuss some interesting contrast of 7 and Z production that can serve as an 
additional marker of the D-brane model. Ignoring the Z-mass (i.e., keeping only transverse 
Z's), and assuming that cross sections x branching into lepton pairs are large enough for 
complete reconstruction to pp — * Z + jet, the quiver contribution to the signal is suppressed 
relative to the photon signal by a factor of tan 2 9w = 0.29. The SM ratio [Z background)/( 
7 background) is roughly 0.92 for processes involving u (or u) quarks, and 4.7 for processes 
involving d (or d) quark. Thus, even if d quark processes are ignored, one obtains a signal- 
to-noise ratio (S/N)z = 0.29/V0.92 = 0.30 (S/N 7 ). Keeping the d quarks will only lead to 
more suppression of (S/N)^ |29(]. This implies that if the high-fc^ photons, as predicted by 
the TeV string model, are discovered at 5a, they will not be accompanied by any significant 
deviation of pp — > Z + jet from SM predictions. This differs radically from the evaporation 
of black holes produced at the LHC. In such a case, production of high-fc^ Z and 7 are 
comparable. The suppression of high-fc^ Z production, whose origin lies in the particular 
structure of the quiver model, will hold true for all the low-lying levels of the string. 
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Appendix I 

The SM background for processes with a single photon in the final state originates in the 
parton tree level processes gq — > jq, gq — > jq and qq — > jg, 



„. da 

2E'- 



d 3 k' 



pp-*yX ijk 



/rlrr 
dx a dx b fi(Xai Q) fj(x b , Q) 2E '-p7J 



(27) 

ij—*"fk 



where x a and x b are the fraction of momenta of the parent hadrons carried by the partons 
which collide, k' (E') is the photon momentum (energy), da/d^k'^^^ is the cross section 
for scattering of partons of type i and j according to elementary QCD diagrams, fi(x a ,Q) 
and fj{xb,Q) are parton distribution functions, Q is the momentum transfer, and the sum 
is over the parton species: g,q = u, d, s, c, b. In what follows we focus on gq —>■ jq, which 
results in the dominant contribution to the total cross section. Corrections from the other 
two processes can be computed in a similar fashion. The hard parton-level cross section 
reads, 



. da 
2E' 



d 3 k> 



gq^-yq 



i^L 1 S [(k + p- k'f] - Y \M\ 2 = j\,^5{2 V . q + q^Y\M\\ 
(2vr) 6 2s [V y ;| 4^' ' (2tt) 2 2s V ^ H y; 4^' '' 

(28) 

where k and p are the momenta of the incoming partons, q = k — k', 's — x a XbS, and 

-q 2 = -t = Q 2 . Here, 



lEM 



A. ry ry rt I O O \~ L 



ijW ° +1^1 , (29) 



where g and e are the QCD and electromagnetic coupling constants, and e q is the fractional 
electric charge of species q. For completeness we note that for qq — > gj, 

'— |M|^VeVJ-^-^l. (30) 




4^ r " "9 yo °^ s + t 

Equation (128]) can be most conveniently integrated in terms of the rapidity y and transverse 
momentum k± of the final photon 

d 3 k' 1 

7T^7 = -d 2 k±dy = Trk±dk ± dy. (31) 
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Considering that the incoming momentum of the gluon is k = x a P\ and that of the quark 
is p = XbP2, we can re- write the argument of the delta function as 

2p . q + q 2 = 2x b P 2 . {x a P x - k') +t = x a x b s - 2x b P 2 . k' + t , (32) 

where Pi and P 2 are the initial momenta of the parent protons. Introducing, k' = k± coshy, 
fcjj = k± sinhy, P x = {■s/s/2, 0, 0, y/s/2), and P 2 = (yfs/2, 0, 0, -y/s/2) we obtain 

P 2 . k' = ^k ± (cosh y + sinh y) = ^ k ± e y (33) 

and 

t = -2k . k' = -2x a \ k± e- y = -y/s~k ± e~ y x a , (34) 



so that 



5(x a x b s — ^fsx b k± e y — y/sx a k± e y ) = - 5(x a x b — x b x±e y — x a x±e y ) 

s 

x a x ± e~ y 



S(x b - -^ , (35) 

» [x a — x± e»j 

r bound x b > implies x a > x± 
leads to a stronger constraint 



s [x a — x± e y ] \ x a — x± e y / 

where x± = k±/y/s. The lower bound x b > implies x a > x± e y . The upper bound x b < 1 



x±e y 



1 — x±e y ' 



(36) 



which requires x±e y < 1 — x±e~ v , yielding x± < (2coshy) _1 . Of course there is another 
completely symmetric term, in which g comes from P 2 and q comes from Pi. Putting all 
this together, the total contribution from gq — > jq reads 

<l^x = zJ2j^I dXa I dXb f '^ Q) &*» Q) J^f s [x a - x ±e y] 



v 1 ,( x a x ± e~ y \ e 2 g 2 e 2 q s + t s \ 

x Ys 6 [ Xb -x-^x^) -3" [— + ¥Tt) • (37) 

With the change of variables z = e y Eq. (|3"7|) can be re-written as 
°£??x = 2 £ / * k± f ±dZ I dx a J dx b f 9 (x a , Q) f q (x b , Q) {2ir?2XaX ^ {Xa _ x±z) 



y * ( ■>-,, - - ^i^lL) -Xi | — — i , m) 



x a x±z~ l \ e g e q (s + t s 



Now, since 



x a -x ± zj 6 \ s s + t 



t y/sk ± e y x ± x ± z _ 



Eq. (|38|) becomes 

; 2 2 /-1/2 /-2 m ax /-l 

<t/;;. ;;. = — — / dx ± / dz dx a f g {x a ,Q) 



pp-* 7 x - U7TS 



*^_l_min " ^min ** <& 



a .trt /ft 



E««/.(a6.o) 



1 fx±z x a \ 

% 2 a \ X a X ± Zj 
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where the integration limits, 



2^raax 
min 



± 



x± 



and 



J : f: 



x±z 



1 — X±Z' 



-1 ' 



(41) 



are obtained from Eq. (l3T)j) . In Fig. [2] we show the QCD background cross section vs k±^ min , 
as obtained through numerical integration of Eq. (|40l) . To accommodate the minimal accep- 
tance cuts on final state photons from the CMS and ATLAS proposals [30], an additional 
kinematic cut, \y\ < 2.4, has been included in the calculation. 



Appendix II 



For the considerations in the present work, the resonant cross section can be safely ap- 
proximated by single poles in the Narrow- Width Approximation, 



r^/Tr 



7T 



7T 



(?- So ) 2 + aV^) 2 iV^ r^ 



5(s - So) 



(42) 



where so = M|. The scattering proceeds through J = and J = 2 angular momentum 
states, with the Mf term in Eq. ( 1201) originating from J = 0, and the t 4 + u 4 piece reflecting 
J = 2 activity. The widths of these two resonances are different, with r J=0 = (3/4) a s M s , 
and r J=2 = (9/20) a s M s [5|]. The average string amplitude square in Eq. (1201) then becomes 



\M(gg^ 91 )\ 2 « Ag*Q'C{N) 



7T 
^5/2 



+ 1 



pj=o 



+ 



t 4 + (t + Sp) 
TV=2 



5(s - s ) 



4g*<?C{N) -?s {14 + f $ 4 + (t + s ) 4 }} S(s-s ). (43) 



<y s Sq 

Thus, the total cross section for single photon production in gluon fusion is given by 
d 3 k' f . f . . , „, . , „, 1 1 



a 



99^19 

pp^^/X 



dx a / dx b fg(x a ,Q)f g (x b ,Q) 



x 



2E i J —J ~~"^v~" t*,'»v~w, Q^y 2?s 

4g 4 Q 2 C(N)— 3 {§s 4 + f p + (?+ s ) 4 ]} 5(s - s ) 



a., s 



5(a; a x fe -x b x_LZ -x a x_L2; x ) 

(44) 



s -^o 



We set Q = M s , which is appropriate for the dual picture of string theory. We are aware that 
for Q ~ M s , the parton distribution functions will receive significant corrections from the 
rapid increase of degrees of freedom. Fortunately, as noted elsewhere [3JJ, at parton center- 
of-mass energies corresponding to low-lying string excitations the resonant cross section is 
largely insensitive to the details of the choice of Q. Plugging r = s Q /s into Eq. fj44l) . we 
obtain 



a 



99^19 
pp—t'fX 



ir k\ dk\ dz 



dx a / dx b f g (x a ,Q)f g (x b ,Q) 



\ cl "^b h _l_ ft _L 

8 n 2 x\x b 



x Ag 4 Q 2 C(N) 



x 5 [ x b 

X f, 



71 S 



{|r 4 + f [(xax^- 1 ) 4 + (-x a x ± z- 1 + r ) 4 ]} 



(45) 
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FIG. 6: Relative contributions of initial state partons (ij = gg, gq, gq, and qq) to 
L /t(»o, Q) fj(ro/x a , Q) dx a /x a , with varying string scale. 



which after integration over xj, leads to 



a, 



pp—>-/X 



g 4 Q 2 C(N) fx ± dx ± dz 



2a s r 4 s 



z 



dx a f g {x a , Q) fg(r /x a , Q) — 

Xa 



X 6 [To 



TqX±Z X a Xj_ 



X a Z 

We now make use of the delta function scaling property. 

TqXj_Z X±X 



) {i^ + f[(x a x ± z- 1 ) 4 + (-x a x ± z- 1 +T ) 4 ]} .(46) 



6 To- 



J (I 



z 



<*(/(*)) 



\fh\ li '-' +H \fkT\ l{z - z - h (47) 



where z± are the solutions to f(z) = 0, 



z± 



X a 

2x7 



i±Wi 



Ax\ 

TO 



(48) 
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Besides 



1 



*fc|/'(*fc)l 
and 



r x ± z± x a x± 



X a Z± 



(49) 



XaX± ? irpjl-^l; (50) 



2± 2 \ V r 

therefore, 

1 = 1 f51) 

*t|/'(*fc)| r ov /l-4xi/r ' l ' 

A straightforward calculation shows that 

f To ( 5 A - 10 xi r + 4 r 2 ) = { I r 4 + f [(x.x^^ 1 ) 4 + {-x a x±z+ 1 + r ) 4 ] } 

+ {|r 4 + f [(xaX^Z 1 ) 4 + {-XaXs.zZ 1 + r ) 4 ]} , (52) 
and hence integration over the z variable yields 
„ 9 _ 8 g*Q 2 C(N) fW ^ 



a 



/ dx ± - l 2 = (5xl-104r + 4r 2 ) 

J x i mi - \/l — 4x i /Tn 



/ — ~ fg(x a , Q) fg(T /x a , Q) , (53) 



X 

'TO 



where the integration range has been derived from the conditions < x\, = Tr,/x a < 1 and 
Ax 2 ± < To, which imply r < x a < 1 and a;j_ jm i n < x± < ^Jtq/2. Finally, integration over x± 
leads to 



f 1 dx 

I — fg(Xa,Q) fg(r /x a ,Q). (54) 

J Tn X a 



99~>19 I qQ C(N) ^ X ±,min ( K 2 a 2 , 4 \ 

G pp^ 1 X - g agT 2 s V L ~ 1° T ~ D T X _L,min + ^ X _L,minJ 

r(XX a 
ro ^a 

Note that all stringy corrections to the pure bosonic cross section given by Eq. (l5~4"j) have 
similar factorizations. An illustration of the relative partonic luminosities of the different 
processes is shown in Fig. [6j 

Appendix III 

We follow the same conventions and notation given in Appendix I for two-body processes 
leading to final states consisting of 7 + jet, with equal and opposite transverse momenta k± 
and p±, respectively. The distribution of invariant masses M 2 = (k' + p') 2 is given by 

da (2vr) 4 f d 3 k' f d 3 p' 



ijk 

x 6(p 2 - M 2 ) ^Y,\ M \ 2 > ( 55 ) 



spins 
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where 

J2\M\ 2 = \M(ij -► 7 A;)| 2 = 64tt 2 ?^ = 16tts 2 ^ , (56) 

spins sj— »T«. 

p 2 = g = (fc' + p ') 2 = 2A;'. p' = 2^;^ - fcjpj + p^ and 

<5 4 (p - A4 - P±) = *(£ - #i - E 2 ) 8{p\\ - k\ - p'{]8(kx_ + p ± ) . (57) 

The integration over d 3 k' d 3 p' can be conveniently re-written in terms of rapidities y\ and yi 
(of the 7 and the jet) and their common transverse momentum, 



d 3 p n 
~2E ~ 2 



dp\ dy , (58) 



where y = \(y\ — 2/2)- Since £?{ = p±coshyi, fen = p_i_sinhyi, I*^ = p_i_cosh?/ 2 , and pj, = 
p_i_ sinh y 2 , a straightforward calculation leads to E[E' 2 — k'np'u = p\ cosh(yi— y 2 ) = p\ cosh 2y. 
Now, using the identity of hyperbolic functions, 1 + cosh 2y = 2 cosh y, we define 

s _ M 2 _ 4pi _ u2 



r = - = = cosh y (59) 

s s s 



so that 



1 / M 2 \ 

5(s-M 2 ) = 5(Aplcosh 2 y-M 2 ) = —— 2 - sU-—-^-) . (60) 

4 cosh y \ 4 cosh yj 

Using 

J d 2 k ± d 2 p ± 5(k ± +p ± ) 5(p 2 ± -M 2 /4cosh 2 y) = n ! dp\ 5(p\ - M 2 /4cosh 2 y) = vr , (61) 
Eq. (|55|) becomes 

1 (8ttM 2 ) J dy x f dy 2 J^ I "<&« /<&& /*(&«, ^) /i(a*,M) — - 



rfM 2 (2vr) 2 4 v " /,j " v yi y y ^w y "w; J3 w«j 4cosh2 



x 5(E-E' 1 -E' 2 )5(p\\-k[-p'{) -* 

dt 



(62) 

ij—ryk 



We now define a — E — E\ — E 2 and 6 = py — A;|, — pj, to perform the change of variables 
A = a + b and B = a — b, such that 6(a) 5(b) = N5(A) 5(B), with normalization iV given by 

/ da db 5(a) 5(b) = f dAdB 3j^L N 5(A) 5(B) = y = 1 ■ (63) 

The new variables can then be explicitly written as {^ } = E ±p\\ — (E\ ± kl) — (E2 ± p|, ) , 
where £±py = {^}, #i ± fcj = P±e ±?/1 = p ± e ±(y+?/) , and £ 2 ±p[ = p ± e ±i/2 = p ± e ±(y ^ } , 

with F = \(yi + 2/2) ■ Putting all this together, the product of delta functions in Eq. (|62|) 
becomes 

#(£ — Ei — £^ 2 ) 5(p\\ — fcj — pj|) = 25(\fsx a — 2p±e Y coshy) 5(\^sx b — 2p±e~ Y cosh?/) 

= 25(y/sx a - Me Y ) 8(Vsx b - Me~ Y ) , (64) 
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and hence integration over the fraction of momenta is straightforward, yielding 



m -\ MT fa*» ^r y £/.(v^) a(^-,m) | 



(65) 

ij—fyk 



Now, if we constrain the rapidities to the interval 2.4 < yi, y% < 2.4 we obtain the invariant 
mass spectrum given in Eq. (124"I) . In addition, note that x a ,Xb < 1, implying — ln(l/y / r) < 
Y < ln(l/y/r). Besides, the cut on the transverse momentum leads to fc_i_, m in < M/2coshy. 
Finally, the Jacobian reads 

dyi dy 2 = *' \ dYdy = 2dYdy, (66) 

d{Y, y) 

and the region of integration is defined by \y\\ = \y + Y \ < 2.4 and 1 2/2 1 = \y — Y\ < 2.4. 
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